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Abstract 



We study the correlation energy associated with the pair fluctuations in 
BCS theory. We use a schematic two-level pairing model and discuss the 
behavior of the correlation energy across shell closures, including the even- 
odd differences. It is shown that the random phase approximation (RPA) to 
the ground state energy reproduces the exact solutions quite well both in the 
normal fluid and in the superfluid phases. In marked contrast, other methods 
to improve the BCS energies, such as Lipkin-Nogami method, only work when 
the strength of a pairing interaction is large. We thus conclude that the RPA 
approach is better for a systematic theory of nuclear binding energies. 
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I. INTRODUCTION 



An important goal of nuclear theory is to predict nuclear binding energies. The mean-field 
approximation certainly provides a good starting point, but the correlation energy associated 
with nearly degenerate configurations can be of the order of several MeV. Since one would like 
an order of magnitude better accuracy for global modeling, the correlations must be treated 
with some care. The correlations associated with broken mean-field symmetries, namely 
center-of-mass localization, deformations, and pairing, are especially important and need to 
be singled out for special attention. The study of such correlation effects has attracted much 
interest, see Refs. |]-^| for recent citations. Popular computational methods [] include the 
generator coordinate method (GCM) 0], variation- with-projection methods (e.g., VAP) ||, 
and the RPA §J. 

In previous work [ID|, we have advocated the RPA approach for deformations and for 



center-of-mass localization. The correlation energy is calculated from the RPA formula 



E, 



corr 



J2^i-Tr(A)\, (1) 



where uii is the frequency of the RPA phonon and A is the A matrix in the RPA equations. 
This method requires solving the RPA equations for each nucleus, which is computationally 
easy if the interaction is assumed to be separable |TT|-[T3|. On other other hand, methods 
such as VAP and GCM are rather complicated to apply. One should also remember that a 
global theory requires a method that is not only systematic but also preserves the relative 
computational simplicity of the mean field approximation. In this respect, the Lipkin- 
Nogami method is a popular candidate for a computationally easy method to go beyond the 
mean field pairing theory |T^-|TBJ 



In this paper, we continue our study of the accuracy of the RPA correlation formula (HI), 
going on to pairing correlations. There is a long history of the RPA treatment of pairing 
correlations. An early study by Bang and Krumlinde |T7j showed that the RPA formula 



reproduces the exact correlation energy rather well in a schematic model. Kyotoku et al. 
|18| compared the leading 1/N behavior of different methods including the Lipkin-Nogami 



method. They found that only RPA gave the exact coefficient in the condensed phase. The 
RPA method has in fact been used in realistic models of deformed nuclei [ IDfl . The RPA 



correlation in the normal phase was studied in Ref. |2(| using the self-consistent version of 
RPA. But perhaps surprisingly in view of the large literature, we have not found any studies 
that specifically compare the RPA with the computationally attractive alternative methods, 
testing the behavior across shell closures and for odd N systems. It is important that the 
method should not introduce spurious discontinuities when the mean field solution changes 



character; otherwise separation energies could not be reliably calculated [21 



One may also ask whether simplified exact solutions may be of use for realistic applications. 
The Richardson solution [^J5| for a pairing Hamiltonian is useful in this respect, but it requires 
a state-independent pairing interaction, which may not be realistic enough. That solution has 
recently been generalized to separable pairing interactions ||, but it is not yet clear whether the 
computational simplicity is preserved. 
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To investigate the applicability of Eq. ([j]) and study the effects of symmetry break- 
ing restoration, in this paper we employ a well-known schematic two-level pairing model 
|i6| , |20|j22| -pT| , and compare its exact solutions to the approximations that are computation- 



ally attractive. The details of the RPA correlation energy are given in the next section. 
They include both the quasi-particle RPA extension of the BCS theory and the RPA for 
the pairing Hamiltonian when the strength is too weak for the mean field approximation to 
support the BCS solution (the "pairing- vibration" regime). The specific application of the 
various methods to the two-level model is given in the following section III. 

II. PAIRING HAMILTONIAN AND THE RPA 

The pairing Hamiltonian and the BCS solution are well-known and we just summarize the 
equations to confirm the standard notation. In this paper we consider a Hamiltonian with an 
arbitrary single-particle term but a pairing interaction whose strength is state- independent, 

# = E^-g£4v, (2) 

j id' 

where €j is a single-particle energy. The number operator Nj for the shell j is given by 

Nj = ( a ]m a jm + a]fh a jm) • (3) 

m>0 

Here, is the annihilation operator for the time-reversal state and is given by a jm 

a ; 



\j-m 



j-m 



. The pairing operator A\ in the Hamiltonian ([J) is given by 



A 1 = V n 1 (4.) 

-rLj / j LL jm u 'jrm V / 

m>0 

and Aj is the Hermitian conjugate of Aj. The exact solutions of the Hamiltonian were 
obtained a long time ago by Richardson and Sherman |4]]. In Appendix A, we solve the 
Richardson equation for a simple, but non-trivial case where there are two pairs in a single 
j-shell. 

The BCS theory is the mean field solution to Eq. (0). We remind the reader that it is 
derived variationally from a ground state wave function of the form 

\BCS)= [] ( u j +Vja) m a] fh )\), (5) 

j,m>0 

where Vj is the usual pair occupation amplitude and Uj satisfies + = 1. This wave 
function is appropriate for even-A r systems, N being the number of particles in a system. 
The extension to odd-iV systems is given at the end of this subsection. Important quantities 
in the variational solution are the chemical potential A and the pairing gap A = G J2j ^j u j v j-> 
Qj = (2j + l)/2 being the pair degeneracy of the j-shell. When the pairing strength G is 
small, the variational equations only have the trivial solution, Vj,Uj = or 1 and A = 0. 
The ground state wave function is thus given by \HF) = I\j. m a ]m\) an d the ground state 



energy is obtained as 
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E HF = {HF\H\HF)= Y fypej-G). (6) 

j -.occupied 

When the strength of the pairing interaction G is larger than some critical value G cr u, the 
variational equations have a non-trivial solution A ^ and the ground state energy is given 
by 

E BCS = (BCS\H\BCS) =2£fi^ - ^ - <?£%4 (7) 

3 G 3 

For odd N systems, one of the particles in a system does not form a pair and blocks a 
level. When the k-th level is blocked, the HF energy, the BCS energy, and the pairing gap 
are modified to 

E HF = Y &i(2e j -G)+e k , (8) 

jioccupied 

E BCS = 2 Y fWj - ^ - G Y Vjv* + e h , (9) 

3 3 

A = GY&jUjVj, (10) 

3 

respectively. Here, Qj = Qj for j ^ k, and Q k — &>k — ^N,odd- The chemical potential A is 
determined so that 



Y^ + 1 = N, (11) 

3 

is satisfied. This formalism is referred to as the blocked-BCS theory. 

A. Random phase approximation 

We next introduce the random phase approximation to compute the correlation energy. 
Let us first consider the RPA in the superfluid phase (QRPA). We refer to Refs. |]7|, |22|j28[ 
for the formulation. The result is the well-known RPA matrix equation 

(-1^) (?)—(?). M> 

where the matrices A and B are given by 



Aij = 2E i 8 itj - G \Jnj(uM + vM), (13) 



Bij = G\IQi ^jQj [u\v] + v 2 iU )) , ( 14) 



respectively. Here, Qj is defined in the previous subsection, and Ej = J (e,- — A — Gvf) 2 + A 2 
is the quasi-particle energy. The RPA excitation operator is given in terms of quasi- 
particle operators 
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^jm ^j^jm Vj(lj mi (l^) 
^jfh = U j a lm v j a jmi (16) 

as 



Q = J2 \ X i X) a Jm«jm - Y j Yl a j™ a j™ / V flj- (17) 



The QRPA correlation energy is given by Eq. ([I]) with the A matrix given by Eq. (|T3|). 
Note that the second term on the right hand side of Eq. (HD is not J2i 2-Ei but Tr(A). This 
point was overlooked in the integral approaches to the correlation energy in Refs. Hlll,[l2|. 

We next consider the RPA in a normal-fluid phase (pp-RPA), which describes pairing 
vibrations. The equation for the pp-RPA can be obtained from the QRPA equation by 
setting A = and — u p — 1, where p and h denote particle and hole states, respectively. 
As for the chemical potential A, there is no definite value for it in the normal fluid phase 
since it can be anywhere between the highest occupied and the lowest unoccupied levels. 
Notice, however, that the pairing vibration describes the ground state in the N ± 2 systems 
and the role of the chemical potential is therefore just to shift the energies by an amount 
±2A. After removing these trivial energy shifts, one obtains 



App* — Idpdp^pi — G \j Vtp\j Vtpi ^ (18) 
A hh , = -2{e h - G)8 h>h , - G^KJK,, (19) 



B p h = Gy Q p \f Q h , (20) 
A p h = B pp > = Bhw = 0. (21) 

This 2(N P + Nh) x 2(N P + Nh) dimensional matrix equation, where N p is the number of 
unoccupied shells and is that of occupied shells, can be decoupled into two (N p + Nh) x 
(N p + N h ) matrix equations as 





and 




(22) 



hu r [ vr \. (23) 




Here, the matrix C is defined as Chh' = A hh i,C p h = 0, and RPA amplitudes are given by 
Xp = X p , = —Yh, X r h = —Xh, and Y p r = Y p . The first equation ( p2f) describes the ground 
state of the iV + 2 system and is referred to as the addition mode, while the second equation 
(f23]) describes the N — 2 system and is called the removal mode. Noticing that Eq. (P is 
a half of a sum of the difference between the RPA and the Tamm-Dancoff approximation 
(TDA) frequencies for each mode of excitation 0, we find the correlation energy for the 
addition mode to be 



C = 5 £^-7V(A) , (24) 
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while that for the removal mode is 

Krr = I (ll^ri-TriC)^. (25) 

The total correlation energy is the sum of these, E corr = E® orr + E r corr . 

Typically, one finds that an RPA mode goes to zero frequency at the mean-field phase 
transition. This is not the case for Eqs. ( p2|) and (p3f), thus giving a discontinuity in the 
frequencies at the transition. However, as we have mentioned above, this is an artifact of 
an awkward choice of the chemical potential A, and it will be the case that the sum of an 
addition and removal mode goes to zero. 



B. Lipkin-Nogami method 



An alternative way to restore the broken gauge symmetry of the BCS approximation is 
to carry out the number projection of the BCS wave function. The Lipkin-Nogami method 



provides an approximate way for number projection. It was first invented by Lipkin [14], and 
was developed by Nogami and his collaborators |15|Jl6|l . Because of its relative simplicity, it 
has been widely applied p^j31| -|34 . 

In the Lipkin-Nogami method, the expectation value (H — XN — \2N 2 )bcs is varied. 
The resultant equations to be solved are given by [[LJJ 

2 ^ a 



G 



E 



li + A* 



2 ^v 2 + 5 NMd = N, 



A 2 = - 



G \ fe tiiufa) fe tUuiV?) - Ei a 



4 4 



fen 



l u 2 i v 2 i 



(26) 
(27) 

(28) 



where e^v^ and Ui are defined as 







-(4A 










-i( 












= 1 - 





(29) 
(30) 

(31) 



respectively. We use the blocked-Lipkin-Nogami prescription for odd iV systems. One of the 
characteristic features of the Lipkin-Nogami method is that the pairing gap A has a finite 
value even in the weak G limit where A is zero in the BCS approximation. The ground state 
energy in the Lipkin-Nogami method is given by 



A 



E L n = 2 ^ fye 3 -vj - - G ^ fyv* - 4A 2 ^j^v 2 + e k 5 N , odd . 



(32) 
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III. TWO-LEVEL PAIRING MODEL 



We now apply the above equations to a schematic two-level model. This model was first 
introduced in Ref. f22|, and has been used in the literature to test several approximations 
16| , pr8|^0| , |23| ^p7[1 . We label the lower and the higher levels 1 and 2, taking t\ = — e/2, e 2 = e/2. 



The Hamiltonian can be numerically diagonalised using the quasi-spin formalism 0. The 
basis states are denoted by (.SiSio; S2S20), where Si and S^ are defined by Si = (Qi — z^)/2 
and Sio = (iVj — f)j)/2, respectively. The latter takes a value of — Si, —Si + 1, • • • , Si. Vi is 
the seniority quantum number. For the ground state of even-iV systems, the total seniority 
Vi + v 2 is 0, while it is 1 for odd-iV systems. The matrix elements of the Hamiltonian read 



(S[S' W ; S 2 S' 20 \H\SiSx ; S 2 S 20 ) — $Si,s' 1 8s 2 ,s 2 je (s 20 — £10 — $s w ,S' 10 $, 

-G (Sx(St + 1) - S 10 (S W + 1) + S 2 (S 2 + 1) - S 20 (S 20 + 1 )) tisio,S',Jsaa,S; 



-GJiS^Si + 1) - S 10 (S 1Q - l)JS 2 (S 2 + 1) - S2 (S 20 + l)6 Sl0 ,s> a +iS. 



S' w +l°S 2 o,S' 20 -l 

-G\J {Sx{S\ + 1) - 5io(5'io + l)\]S 2 {S 2 + 1) - S 20 (S 20 - 1)8s 10) s' 10 -iSs2o,s' 20 +i^ ■ (33) 

We first consider a symmetric two-level problem i.e., Qi — Q 2 — and assume the 
number of fermionic particle is N = 2Q. Equations for several approximations can be solved 
analytically for such a system. The gap equation in the BCS theory leads to a pairing gap 
of 



together with 



A = WG»n*- T , (34) 



*-H( 1 + 2gnH' (35) 

u\ -v^Ul-J-)^ u\ (36) 



2 V 2GQ 

X = — - Gv\ (37) 



where e is defined as e = 2Qe/ (2Q — 1). 

From Eq. ([34]) , one can find the critical strength of the phase transition to be G CT n = 
e/ (2Q — 1). For G larger than G cr i t , the system is in the superfluid phase, and the matrices 
A and B for the QRPA equation are given by (see Eqs. (|H|) and ([14])) 

A2 * (38) 



An 


= A 22 


= Gn 

2 


A l2 


= A 21 


A 2 

2GQ' 


B n 


= B 2 2 


A 2 
2GQ 2 


B\2 


= B 2 i 


= Gtt 





GO? 2GQ' 

(39) 

+ 2^ ^ 
A 2 

2«r (41) 
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The solutions of the QRPA equation are Too = and y4A 2 — 2A 2 /f2, from which the 
correlation energy is computed as 

A 2 A 5 



+ — . (42) 

When the strength of the pairing interaction G is smaller than G cr u, the system is in the 
normal fluid phase. The mean field energy is then evaluated as Ehf = —^Q — GQ. For the 
present two-level model, the A, B, and C matrices for the pp-RPA are just numbers and are 
given by A = e — GQ, B = GQ, and C = e — Gfl + 2G, respectively. The frequencies for the 
addition and the removal modes are then found to be 

hu a — —G + VTTGVe + G- 2GVL, (43) 
%uj r = G + Ve + GVe + G- 2GVL, (44) 

respectively. The total correlation energy is thus obtained as 



E corr = Ve + GVe + G- 2Gtt - (e — GQ + G). (45) 

Figure 1 shows the RPA frequencies for each modes of excitation as a function of G, for 
Q — 8. As we noted in the previous section, we see that because of the chemical potential 
the pp-RPA frequencies do not match with the QRPA frequencies at the critical point of 
phase transition from a normal-fluid to a superfluid phases. Figure 2 compares the ground 
state energy obtained by several methods. The solid line is the exact solution obtained by 
numerically diagonalizing the Hamiltonian. The dashed line is the ground state energy in 
the mean field BCS approximation. It considerably deviates from the exact solution through 
the entire range of G shown in the figure. The dot-dashed line takes into account the RPA 
correlation energy in addition to the mean field energy. It reproduces very well the exact 
solutions, except in the vicinity of the critical point of the phase transition. Around the 
critical point, one would need to compute the correlation energy with some care, using e.g., 
the self-consistent RPA discussed in Ref. |2(| which removes the cusp behaviour of the RPA 



frequency around the critical strength. 

It is interesting to compare the present mean field plus RPA approach with the Lipkin- 
Nogami method. The equations for the Lipkin-Nogami method for the two-level model were 



solved in Ref. [HJ. The results are given by 

A 2 = (Gn) 2 (l - k 2 ), (46) 
v\ = u\ = (1 + «)/2, (47) 
u \ = v l = (1 _ «)/ 2 , (48) 

4A --^ (2,-l)(f-^) ' ^ 
where k is the physical solution, which satisfies < k < 1, of an equation 

2(1 - Q)k 3 + (2fi - l)nk 2 + (2Q - l)k - (2Q - 1)k = 0, (50) 

and k is defined as e/2GQ. The ground state energy in the Lipkin-Nogami method given by 
Eq. (|3~2"1) is denoted by the thin solid line in Fig. 2. Although it reproduces the exact results 
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for large values of G, it deviates significantly from them for small values. This behaviour is 
consistent with the numerical observation in Ref. [^4| as well as the result of Ref. [18| where 
it was found that the Lipkin-Nogami method is only correct in the limit of a strong pairing 
force. In marked contrast, the ground state energy in the BCS plus QRPA coincides with 
the exact solution at the leading order of an expansion in 1/Q for any strength of the pairing 
interaction as long as it is larger than the critical value [|l8|j . 

From Fig. 2, it is unclear whether the Lipkin-Nogami method or RPA are more accurate 
for purposes of a global theory of binding. So we now consider a more realistic situation, 
varying the particle number N rather than the interaction strength G. We consider the 
paring energy in oxygen isotopes, taking the neutron lp and 2s- Id shells as the lower and 
higher levels of the two-level model. The pair degeneracy Q thus reads Q\ = 3 and ^2 = 6, 
and the number of particle in a system is given by iV = A - 8 - 2 for the A nucleus. We 
assume that the energy difference between the two levels e is given by e = AlA^ 1 ^ and the 
pairing strength G = 23 /A. The upper panel of Fig. 3 shows the ground state energy as a 
function of A. In order to match with the experimental data for the 16 nucleus, we have 
added a constant —72.8 MeV to the Hamiltonian for all the isotopes. The exact solutions 
are denoted by the filled circles. The deviation from the BCS approximation (the dashed 
line) is around 2 MeV for even A systems and it is around 1.2 MeV for odd A systems. 
This value varies within about 0.5 MeV along the isotopes and shows relatively strong A 
dependence. One can notice that the RPA approach (the dot-dashed line) reproduces quite 
well the exact solutions. On the contrary, the Lipkin-Nogami approach (the thin solid line) 
is much less satisfactory and shows a different A dependence from the exact results. The 
pairing gaps A in the BCS approximation and in the Lipkin-Nogami method are shown 
separately in the lower panel of Fig. 3. For the Lipkin-Nogami method, we show A + A2, 
which is to be compared with experimental data |TJ|,|TB[. The closed shell nucleus 16 and its 
neighbour nuclei 15>17 have a zero pairing gap in the BCS approximation, and the Lipkin- 
Nogami method does not work well for these nuclei, as can be casted in Fig. 2. On the 
other hand, the RPA approach reproduces the correct A dependence of the binding energy. 
Evidently, the RPA formula provides a better method to compute correlation energies than 
the Lipkin-Nogami method, especially for shell closures. 



IV. SUMMARY 

Returning to our initial motivation, we seek a computationally tractable way to include 
pairing effects in a global model of nuclear binding energies, going beyond the BCS theory. 
Two attractive possibilities are the Lipkin-Nogami method and the RPA, in particular if 
the pairing interaction has a separable form. In this paper, we used a solvable two-level 
pairing Hamiltonian to show that the RPA formula for the correlation energy reproduces 
well the exact solutions both in a normal-fluid and a superfluid phases. On the contrary, 
the Lipkin-Nogami method is considerably less accurate for weak pairing, and therefore 
is not suitable in transition regions. As a consequence, the Lipkin-Nogami method fails 
to reproduce the correct mass number dependence of the binding energy around a shell 
closure. The correlation energy associated with the number fluctuation for the neutron 
mode in O isotopes was shown to be order of 2 MeV and has a relatively strong mass 
dependence. Although the correlation energy is small compared with the absolute value 
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of the ground state energy, this suggests that including the correlation energy in the RPA 
provides a promising way to develop a better microscopic systematic theory for nuclear 
binding energies. 
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APPENDIX A: RICHARDSON SOLUTION 

For a system where the number of pairs is N pair , i.e., a 2iV pa j r -fermion system, the ground 
state energy of the Hamiltonian (0) is given by [|,|J 

E gs = £ *A, ( Al ) 

A=l 

where z\ are N pair solutions of N pair coupled equations given by 

Y— ^ = -• (A2) 

y 2 Ci - z x y z y - z x G K J 

The prime in the summation of A' means to take only those A' which are different from A, 
and Qj = (2j + l)/2 is the pair degeneracy of the j-shell. In general, z\ are complex. 

For a two pair system (N pa i r —2) in a single j-shell, setting 6j = and Qj = Q, the 
Richardson equation reads 



Q 2 1 



Z\ Z2 — Z\ G 

Q 2 1 



Z2 Z\ — Z2 G 

The solutions of these equations are found to be 



(A3) 
(A4) 



z = -G(n-l)±t^G 2 (n~l), (A5) 

from which the ground state energy reads E gs = —2G(Q— 1). This result coincides with the 
solution obtained using the seniority scheme, 

E(N pmr ) = -G || (f + l) " |(2AW - fi ) 2 + \(2N pair - ft)} . (A6) 
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0.05 0.1 0.15 

G/e 

FIG. 1. RPA frequencies as a function of the strength of the pairing interaction G for a 
two-level system. The pair degeneracy ft is set to be 8. The strength of the pairing interaction G 
is given in the unit of e. 
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G /£ 

FIG. 2. Comparison of the ground state energy E gs obtained by several methods. The 
solid line is the exact numerical solution, while the ground state energy in the Hartree-Fock (BCS) 
approximation is denoted by the dashed line. The dot-dashed line takes the RPA correlation energy 
into account in addition to the HF (BCS) energy. The results of the Lipkin-Nogami method are 
denoted by the thin solid line. 
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FIG. 3. The ground state energy E gs (the upper panel) and the pairing gap (the lower panel) 
for oxygen isotopes estimated with the two-level model as a function of the mass number. The 
exact reults are denoted by the filled circles, while the meaning of each line is the same as in Fig. 
2. For the pairing gap in the Lipkin-Nogami method, A2 is added to the pairing gap A. 
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